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ABSTRACT 

Power-law distributions are a near universal feature of energetic particle spec- 
tra in the heliosphere. Anomalous Cosmic Rays (ACRs), super- Alfvenic ions in 
the solar wind and the hardest energetic electron spectra in flares all have energy 
fluxes with power-laws that depend on energy E approximately as E~^'^. We 
present a new model of particle acceleration in systems with a bath of merg- 
ing magnetic islands that self-consistently describes the development of velocity- 
space anisotropy parallel and perpendicular to the local magnetic field and in- 
cludes the self-consistent feedback of pressure anisotropy on the merging dynam- 
ics. By including pitch-angle scattering we obtain an equation for the omni- 
directional particle distribution f{v, t) that is solved in closed form to reveal 
(corresponding to an energy flux varying as E~^'^) as a near-universal solution 
as long as the characteristic acceleration time is short compared with the char- 
acteristic loss time. In such a state the total energy in the energetic particles 
reaches parity with the remaining magnetic free energy. More generally, the re- 
sulting transport equation can serve as the basis for calculating the distribution 
of energetic particles resulting from reconnection in large-scale inhomogeneous 
systems. 

Subject headings: acceleration of particles — magnetic reconnection — Sun: 
corona — Sun: flares 



1. INTRODUCTION 

Accelerated particles with power-law spectra are a nearly universal feature of helio- 

spheric plasmas and also characterize the co smic ray spectrum. Anomalous Cosr aic Rays 

(ACRs) Jstone et al.ll2008l : lOecker et al.lboioh . super-Alfvenic ions in the solar wind (iFisk fc Gloecklerl 



^University of Maryland, College Park, MD, USA; drake@umd.edu, swisdak@umd.edu 
^Astronomy Department, Boston University, Boston, MA, USA ; rfermo@bu.edu 



-2- 



20061 ) and the hardest energetic electron spectra in flares (IHolman et al.l 120031 ) all have en- 
ergy fluxes with power- laws that depend on energy E approximately as E~^'^ . An important 
question is whether there is a common acceleration mechanism in these very disparate envi- 
ronments. 

A range of acceleration mechanisms have been proposed to explain the spectra of en- 
ergetic electrons (up to several MeV) and ions (up to several GeV) in imp ulsive flares. 



inclu d ing the reconnection process it s elf and reconnection-driven turbulence (IMiller et al. 



1997; Dmitruk et al. 2004 



Liu et al.l l2006l : IZharkova et al.l 120111 ). Signiflcant challenges 



have been to explain the large numbers of accelerated electrons an d the surprising effi- 



ciency of the conversion of magnetic e nergy to the energetic particles (ILin fc Hudson! Il971 



Emslie et al.ll2005t iKrucker et al.ll2010l ). The single x-line model of reconnection in flares, in 
which electrons are accelerate d by parallel elect ric flelds, can not explain the large num- 
ber of accelerated electrons (IMiller et al.l Il997l ) . O n the other hand both observations 



(ISheelev et al. 



Drake et al. 



2011 



200 



2006b 



Fermo et al. 



a 



Savage et al.ll20 12 ) and modeling ( K 



Onofri et al. 2006; Oka et al. 2010 



lem 



1994 



Shibata fc Tanuma 



Huang et al. 



2011 



2001 



Daughton et al. 



20121 ) suggest that reconnection in flares involves the dynamics of large 



numbers of x-lines and magnetic islands or flux ropes. In magnetohydrodynamic (MHD) sim- 
ulations of multi-island reconn ection test particle s rapidly gain more energy than is available 
in the driving magnetic fleld ( lOnofri et al.l 120061 ). Thus, developing a model of particle ac- 
celeration in a multi-island reconnecting environment with feedback on the driving fields is 
the key to understanding flare-produced energetic particle spectra. 

The seed population of ACRs ar e interstellar pickup particles sin ce their composition 
matches that of interstellar neutrals (ICummings fc Stond Il996l . 120071 ). However, the con- 
ventional idea that they are accelerated at the termination shock (TS)( lPesses et al.lll98ll ) 
was called into question when the Voy agers crossed the ' 
the ACR spectrum did not peak there ( Stone et al. 20051. 



?S ari d found that the intensity of 
2008h. A possible alt e rnate sourc e 



is ma gnetic reconnection of the sectored heliosheath (iLazarian fc Opherl l2009l : iDrake et al. 
20101 ). Simulations of reconnection in the sectored fleld region revealed that t he dominant 



heating mechanisrn was Fermi reflection in co ntracting and merging islands (IDrake et al. 
2010l : iKowal et al.l 1201 ll : ISchoeffler et al.ll201ll ). Because contraction increases the energy 



parallel to the local magnetic fleld and reduces the perpendicular energy, the heating mech- 
anism drives the system to the flrehose stability boundary a = 1 — {/3\\ — (3\)/2 = where 



reconnection is throttled because the magneti c tension drive is absent (IDrake et al.ll2006b 
20 id : lOpher et al.lboill : ISchoeffler et"aDl201l[ ). 



While multi-island simulations have revealed that Fermi reflection in contracting islands 
controls energy gain and drives the system to the marginal flrehose condition, a rigorous 
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model for particle acceleration in such a multi-island system has not yet been developed. 
The Parker eq uation does not describe particle acceleration in nearly incompressible systems 
( iParkerl Il965l) and extens ions do not account for the geometry of reconnection and island 
merging (lEarl et al.lll988[ ). In the present manuscript we explore particle acceleration in a 
bath of merging magnetic islands with a particle distribution function f{v\\,v±) that accounts 
for the velocity space anisotropy along (fy) and across {v±) the local magnetic field and 
includes a phenomenological pitch-angle scattering operator. Thus, the pressure anisotropy 
can be directly evaluated and the feedback on island merging calculated. 



2. Particle dynamics during island merger 



We develop a probabilistic model of particle acceleration in a bath of mergi ng 2-D mag 



netic islands with a distribution of magnetic flux ip and area A given by g{ip, A) (IFermo et al. 



20101). The development of structure in 3-0 may ultimately be important and should be ad- 
dressed but observations (jPhan et al.ll2006l ) and simulations ( iHesse et al.ll200ll ) suggest that 
at the largest scales reconnection is nearly 2-D and this limit is therefore a reasonable start- 
ing point. We first calculate the particle energy gain during the merging of two circular 
islands of radii ri and r2 with rj = Aj/ir as shown in the schematic in Fig. [H Merging 
leads to a sin gle island of area A f = Ai + A2 and with magnetic flux ipf given by the larger 
of ipi and ip2 (IFermo et al.ll2010l ). The reduction of energy by the factor {ipf +ip2)/ip} results 
from the shortening of the field lines as merging proceeds. Thus, energy release takes place 
not at the merging site, but as reconnected field lines contract after merger. As long as the 
kinetic-scale, boundary layer where reconnection occurs is small compared with the island 
radii, the dominant energy exchange with particles takes place on the closed, reconnected 
field lines that release magnetic energy as they contract. 

We take advantage of two adiabatic invariants, the magnetic moment /i = ■mvj_/B and 
the parallel action ^ v\\di, which are constants if the gyration time of particles around the 
local magnetic field and their circulation time around islands are short compared with the 
merging time. The former describes the reduction in v± as B decreases and the latter the 
increase in v\\ as i decreases. The parallel action invariant is valid for velocities that exceed 
the local Alfven speed, which im plies that a seed heating mechanism is n eeded for low (3 
systems such as the solar corona (jPrake et al.l l2009al : iKnizhnik et al.l I2OIII ) but not in high 
(3 systems such as the sectored heliosheath. To calculate i, we first calculate the merging 
velocity r^ep of two islands with differing radii and magnetic fields, Vsep = ri + r2 = —ip{,Bi + 
B2)/{BiB2), since merging magnet ic islands reconnect th eir magnetic flux at the same rate. 
The reconnection rate is given by ICassak &: Shayl (120071 ). ip = 2V12B1B2/ {Bi + B2), with 
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Vi2 = ^r^/c^uB}^2/^^U>J where ~ 0.1 is the normahzed rate of reconnection and = 
1 — 47r(p|| — p±)/ {BiB2). Thus, f-sep = — 2Vi2 and Vu is the island merging velocity. The 
rate of line shortening can now be calculated from the total merging time (ri + r2)/(2Vi2) 
and the difference between the intial field line length as merging starts and the final length 
using area conservation, i = — 27r/ii2Vi2 with hi2 = 2(ri + r2 — \/rf + r^) / (ri +r2). Parallel 
action conservation then yields an equation for f y , 

dv\\ _ hi2Vi2 , . 

To obtain the corresponding equation for v±, we use the conservation of magnetic flux and 
area as a flux tube contracts so that B/^ is constant. Therefore, from /i conservation v']_/^ 
is also constant and 

dvl _ 2 ^12^^12 

dt Ti + r2 

Thus, the perpendicular energy goes down during island merger as the parallel energy in- 
creases. 



3. A kinetic equation for particle acceleration during island merger 

From the energy gain of particles in merging islands we can formulate a model of particle 
acceleration in a very long current layer of length L. Particles are injected into the bath of 
interacting islands in the current layer from upstream as each individual island grows due to 
reconnection of the upstream field. They then undergo acceleration in the merging islands 
until they are convectively lost. The rate of injection of particles is given by the upstream 



distribution function f, 



111 



islands At (iFermo et al. 



v) ti mes the integrated rate of area increase of all of the magnetic 



2010f ). 



/•OO /"OO 

At = 2'KerCAup I / dAdtfj rg{tlj, A), (3) 
Jo Jo 

with the island radius given by r = yJ~Afn. The probability of two islands of radii ri and 
T2 merging is given by their overlap probability 4rir2/L^. Using the conservation of phase 
space volume and summing over all merging islands in the layer, we obtain a differential 
equation for 



where 



D I I ^T^O 4rir2/ll2V^12 

R= I I dld2g,g2jj^^-^. (5) 
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with di = dAidipi. In earlier simulations of multi-island reconnection strong pressure anisotropy 
with p \\ > jO i within the core of mergin g islands was limited by anisotropy instabilities 



( iDrake et al.l l2010l : ISchoeffler et al.l |2011| ) so we have included a phenomenological pitch- 
angle scattering operator of strength u that acts on the angle ( = v\\/v to reduce anisotropy. 
Importantly, the drive R is independent of the particle velocity. It depends on the pressure 
anisotropy through the merging velocity V12. The integral over islands includes only interac- 
tions for which V12 is real (ai2 > 0). To estimate the scaling of R we note that Nt = J diQi 
is the total number of islands in the layer, so for densely packed islands we can define a 
characteristic island radius = L/2Nt. Thus, R ~ erCA/rN- Of course, R can be much 
smaller if ai2 in the expression for V12 approaches zero. 

If / were isotropic and therefore only a function of f , the energy drive operator in Eq. (jl]) 
would vanish when averaged over the angle C,. In this limit t here is zero n et energy gain. 



consistent with Parker's equation in the incompressible limit (jParkerl Il965l ). Equation (jl]) 
is an equi-dimensional equation and therefore has no characteristic velocity scale. Solutions 
therefore take the form of power-laws. An important property of such an equation is that 
the fluid moments of a given order completely decouple from those of differing order and 
their solutions can therefore be readily obtained from Eq. (jl]) in closed form. Specifically an 
equation for p\\ and p_L can be obtained so that ai2 in the energy drive R can be evaluated 
explicitly. Thus, the feedback of energetic particles on the dynamics of reconnection can 
be computed. In the case of no source, sink or scattering, for example, Eq. (jl]) yields 
dp\\/dt = 2Rp\\ and dp±/dt = —Rp± so that p\\ and p± increase and decrease in time, 
respectively, but the total energetic particle pressure p increases, dp/dt = {2R/3){p\\ — p±). 

Instead of directly evaluating the full moments of Eq. (jlj), we simplify the equation by 
ordering the magnitudes of the rates R, v and c^//^. Since the scattering represented by 
V arises from the pressure anisotropy driven by contraction, we argue that u ^ R. On the 
other hand as the spectrum begins to saturate at firehose marginal stability R is reduced 
and z/ increases. We therefore take z/ ^ i? ~ er-CA/rN ^ ca/L, where the latter follows 
because L ^ r^. The large z/ assumption allows us to solve Eq. (jlj) by expanding / in 
a series of Legendre polynomials / = Ylj ■^jiOfji'^) where Pj is the jth order Legendre 
polynomial. By the symmetry in v\\, fi is zero. The equation for /2 follows from balancing 
the reconnection drive acting on /q with the scattering operator acting on /2-P2(C)! f2{v) = 
— {Rv/Qv)dfQ{v)/dv. By averaging Eq. (jlj) over the scattering term vanishes and the 
energy drive term acting on /2-P2(C) yields an equation for /o(f ), 

dfo i?^ 1 (9 4 (9 CAup r \ r (a\ 

This equation is again of equi-dimensional form and has power-law solutions whose indi- 
vidual moments can be calculated. Evaluating the density in steady state, for example, by 
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integrating over velocity, the drive term vanishes and the total number of particles undergo- 
ing acceleration rix is given by rix = Axriup, where = AtL/caup is the integrated area of 
all of the islands in the layer. The firehose parameter needs to be self-consistently evaluated 
and for this we need 



1 /•! /-cx) 1 

P± = ^t^Cy^ dv27:v'mivl - -vl)Uv)P2{(). (7) 



Using the expression for /2 and noting that f y — v'j_/2 = v'^P2{C), we obtain 

R '•^ 



P\\ - P± 



QAtv 

where po is the isotropic pressure calculated from /o 



/•oo 

/ dv A'Kmv'^ fQ{v) = —po, (8) 



^ Pup 

1 - R^L/3cAup^' 



The firehose parameter becomes 



47rp^ ^ no) 

52 1- R^L/3cAupi^' ^ ' 

where B is the average island magnetic field strength based on the sum in Eq. A key 
feature of Eq. ( ITOl) is its singular behavior when 5 = R'^L/ScAupJ^ = 1- This singularity can 
be understood from the power-law solutions to /o, which describe its behavior at energies 
greater than that of the source fup- Taking /o oc v~^, from Eq. ([6]) we obtain 7(7 — 3) = 10/5 
so that when 5 = 1,7 = 5. The second solution, 7 = —2, corresponds to divergent behavior 
and must be rejected. The singularity in Eq. fllOp therefore arises when /q oc and 
corresponds to a divergence of the pressure integral. Thus, it is clear that the requirement 
that the pressure be bounded requires that 7 > 5 or 5 < 1. In deriving Eq. for /o we 
have assumed large scattering so that u ^ R. On the other hand, since island contraction 
drives the anisotropy, we argued previously that u ~ erCA/r^. Thus, 6 ~ erL/rjsi ^ 1. 
Namely, the acceleration rate should always exceed the system convective loss rate since L 
is much larger than the characteristic island size. The resulting divergence of the pressure 
can only be avoided if the reconnection drive R oc a/o is reduced by its approach to firehose 
marginal stability, which forces R <^ u. Unless Pup is very large, the only way that the 
firehose condition in Eq. (fTOl) can be reached is if 5 1 or 7 ^ 5 and /o oc v^^. 

The total energy content Wq = 3po/2 of this high energy tail can be directly calculated 
from the pressure in Eq. using a ~ and 5^1, 



Wo = ^\/3iyL/cAup. (11) 
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Thus, depending on the level of scattering, the total energy density of the energetic particles 
is of the order of, or somewhat greater than, the remaining magnetic energy. In a system 
with low initial /3 equipartition between energetic particles and magnetic field is energetically 
accessible. In a system with high initial /3 equipartition can only be reached if the system is 
open such that energetic particles can access additional sources of magnetic free energy. 



Discussion 



We have derived a general equation (Eq. (jl])) for particle acceleration in a bath of 
merging magnetic islands in a large 1-D current layer. We demonstrated that the E~^'^ 
spectrum is a nearly universal feature of a multi-island reconnecting system for all values 
of initial /3 as long as the nominal acceleration time of energetic particles is shorter than 
their loss rate. This is the correct limit as long as the characteristic magnetic island radius 
is much smaller than the system scale size L. We argue therefore that the widely observed 
E^^'^ spectrum in the heliosphere is a natural consequence of multi-island reconnection. The 
total energy content of this E~'^'^ spectrum reaches parity with the remaining magnetic field 
energy in the system. 

Equation (jl]) can be readily generalized to a 2-D system by replacing the factors 2rj/L 
by Arf/L"^ in the drive term R. The estimate for the scaling of R is unchanged. The model 
loss term CAupf /L should also be replaced by the convective loss rate u ■ V/ with u the 
convective velocity of the system. The arguments leading to the / oc also apply to 
the 2-D equations. In a system in which the driver R is spatially non-uniform the 2-D 
version of Eq. (jl]) could then be numerically solved for the spatial distribution of energetic 
particles from reconnection. The impact of the finite structure of ma gnetic islands that 



might develop in the third direction rema ins an important open issue (jOnofri et al.l 12006 
Schreier et aPboid [Paughton et al.lboilj ). 



There have now been several published sim ulations of particle acceleration and assoc i- 
ated spectra in 2-D multi-current layer systems (jDrake et al.ll2010l : iDrake fc Swisdakll2012l ). 
We can compare the spectra predicted from our equation with the results of those simula- 
tions. Since the simulations were doubly periodic, there was no convective loss. Further, the 
pressure anisotropy was strong so we consider the non-scattering limit of Eq. (jlj) in which 
the source and loss terms are discarded. The exact solution for / is given by 



fivtvi,t) = fivie 



-2G{t) „,2 G(i) 



0), 



(12) 



where G{t) = J^drR{T). This is consistent with exponential growth of the effective parallel 
temperature and an exponential decrease in the perpendicular temperature. The omnidirec- 
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tional distribution function can be computed numerically for any specified initial distribution 
function for comparison with simulation data. The comparison is made with a system with 
sixteen initial current layers in a 409.6(ij x 204.8(ij domain, where di = c/upi is the ion in- 



ertial length (iDrake et al.ll2010l ). In Fig. |2]we show the magnetic field strength at late time 
(t = 100f2~^) in the simulation after islands on adjacent current layers have overlapped. The 
typical island radius rjv at this time is around 15di. The characteristic acceleration rate 
R ~ ^tCa/'^n ~ 0.007f2^^, where Vt^ is the ion cyclotron frequency. Reconnection remains 
strong for a total time of around lQQVt~^ when the pressure anisotropy shuts off reconnection. 
Thus, the integrated acceleration rate is G ~ 0.7. The comparison between the model and 
the simulation data is shown in Fig. |3l The particle energy spectrum from the simulation is 
shown in the initial state and at t = 200fi^^ in the solid lines. Note that the initial state is 
not a simple Maxwellian because of the shift in the ion velocity distribution that is required 
in the current layers. The fit of the initial spectrum with a single Maxwellian, shown in the 
dot-dashed line in Fig. [3] therefore matches the low energy portion of the spectrum very well 
but underestimates the number of particles at high energy. The late time energy spectrum 
from the solution given in Eq. ( !T2|) . after integration over the angle C, is given by the dashed 
line in Fig. |31 The best fit corresponds to G = 0.82 rather than the estimate of 0.7. The 
model reproduces the overall late-time energy spectrum very well but modestly overesti- 
mates the number of particles in the high energy tail. This is probably because the ions 
in the initial spectrum have thermal speeds that are sub-Alfvenic so the Fermi acceleration 
of the low energy io ns is delayed until they gain sufficient energy in reconnection exhausts 
flPrake et allboOQbh . 



Observations in the qui et solar wind have reve aled that the super- Alfvenic ions display 



an fiy) oc distribution (IFisk fc Gloecklen 120061 ). It has been suggested that solar wind 
turbulence would be dissipated in reconnection current layers (IServidio et al.l 120091 ) and 
therefore that reconnection is an important dissipation mechanism in the turbulent solar 
wind. Solar wind observations also reveal that the pressure anisotropy bumps against the 
firehose thresho ld in some regio ns and that there are enhanced magnetic fluctuations at 
these locations flBale et al.ll2009l ). There are therefore mechanisms in solar wind turbulence 
driving anisotropy and the anisotropy is limited by enhanced scattering. Finally, the direct 
observations of reconn ection events in the s olar wind reveal heating but no localized regions 
of energetic particles (iGosling et al.l 120051 ). This is consistent with our picture that the 
energetic particle spectrum is not produced at a single x-line but requires that the ions 
interact with many reconnection sites. 

The spectrum of energetic electrons in impulsive flares are not measured in situ and must 
be inferred from chromospheric x-ray emission. Nevertheless, the energetic pa rticle fluxes do 
occasionally reveal spectra as hard as which corresponds to / oc ( iHolman et al. 
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20031 ). In recent over-the-limb observations of flares in which the reconnection region high 
in the corona can be directly diagnosed, it was found that all of the electrons in the accel- 
eration region became part of the energetic comp onent, indicating th at all electrons in the 
region of energy release underwent acceleration (IKrucker et al.l |2010| ) , which is consistent 
with our model. The /3 of these electrons was of order unity, which is also consistent with 
our predictions. 

Whether the sectored heliosheath magnetic field has reconnected remains an open issue 
because the Voyager magnetome ters are at the limit s of their resolutions at the magnetic 
field strengths in the heliosheath ( iBurlaga et al.ll2006[ ). Large drops in the energetic electron 
and ACR population as Vo yager 2 exited from the sectored zone are consistent with recon- 
nection as the ACR driver f lOpher et al.ll20111). The spectral i i idex of the ACR pa rticle flux 
measured at Voyager 1 is slightly above 1.5 (jstone et al. 2008 : Decker et al. 2010 ). Further, 
the integrated energy density of the measured ACR spectrum between 1 and lOOMeV^ and 
is comparable to that of the magnetic field, which has a magnitude of around 0.15nT. This 
is again consistent with the predictions of our model. 

The equations presented here were derived in the non-relativistic limit. However, the 
ideas can be easily extended to the case where the particles are relativistic but where re- 
connection itself is non-relativistic. We express the distribution of particles in terms of the 
particle momentum p. For the pressure integral to remain bounded 7 > 4 for power-law 
distributions with fo{p) oc p~"'. The resulting particle fiux per unit energy interval F is given 
by F oc p'^fip) (X p'^~"' . Thus, the spectrum of the fiux in the strongly relativistic limit should 
scale as p^^. 



This work has been supported by NSF Grant AGS1202330 and NASA grants APL- 
975268 and NNX08AV87G. 
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Fig. 1. — Schematic of the merger of two islands of differing radii and magnetic field strength. 
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Fig. 3. — The spectra of ions from the simulation of Fig. [2] (sohd hnes) at t = 0, 200^2^/ and 
from the model (dashed and dot-dashed lines). 



